Spectral function of a quarter-filled one-dimensional CDW insulator 
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We consider a one-dimensional charge density wave (CDW) 
insulator formed by Umklapp processes in a quarter-filled 
band. The spectrum of the model consists of gapless, un- 
charged excitations carrying spin ±1/2 (spinons) and gapped, 
spinless excitations carrying charge =Fe/2 (solitons and anti- 
solitons) . We calculate the low-energy behaviour of the single- 
electron Green's function at zero temperature. The spectral 
function exhibits a featureless scattering continuum of two 
solitons and many spinons. The theory predicts that the gap 
observed by Angle Resolved Photoemission (ARPES) is twice 
the activation gap in the dc conductivity. We comment on 
possible applications to PrBa2Cu307 and to the Bechgaard 
salts. 

PACS numbers: 71.10.Pm, 72.80.Sk 



The TMTSF and TMTTF families of quasi one dimen- 
sional (ID) organic conductors have attracted much at- 
tention over the last 30 years. They feature a rich variety 
of 3D ordered ground states, including spin and charge 
density waves, spin-Peierls insulators and superconduc- 
tors. Recently much of the attention has focussed on 
the ID phases at low temperatures, which are be- 
lieved to realize exotic strongly correlated Mott or CDW 
insulating states. With regard to a theoretical descrip- 
tion of these phases there are two main points of contro- 
versy. Firstly, it is not clear how "one-dimensional" in 
particular the TMTSF salts are in the relevant regime of 
energies/temperatures, and secondly it is unknown how 
strong the effective electron-electron interaction at low 
energies is. Persuasive arguments in favour of strong in- 
teractions, giving rise to a small value of the Luttinger 
liquid parameter K c , are given in The authors of 

quote an estimate of K c sa 0.22. At such values of K c 
the single-particle tunneling between the chains should 
be strongly suppressed, leading to the enhancement of ID 
effects. Perhaps, the strongest evidence in favour of such 
a scenario comes from ARPES measurements {|,[|. Ex- 
periments done by different groups not only show the ab- 
sence of quasi-particle peaks, but find no evidence of any 
dispersing features in the single-electron spectral func- 
tion. Quite similar behaviour has recently been observed 
in the quarter-filled Cu — O chain material Prl^CuaCv 
("P123") §§. 

The great difficulty in determining dynamical response 
functions theoretically is the presence of Umklapp scat- 
tering processes, which dynamically generate a spectral 
gap. In this work we use a method based on the inte- 



grability of the low-energy theory to calculate the single- 
particle Green's function for the first time. In the Bech- 
gaard salts there are two separate mechanisms that lead 
to Umklapp scattering: (1) "double" Umklapp processes 
due to the commensurate band filling 1/4 |^||. These 
generate a gap only for strong interactions (K c < 0.25). 
(2) a small dimerization, which halves the Brillouin zone 
and gives rise to "single" Umklapp processes ]l0| , pd| . 
These open a gap already for weak interactions K c < 1 
but their coupling constant is proportional to the dimer- 
ization and thus small. It is an open question which of 
these two mechanisms dominates in the Bechgaard salts. 
Here we develop a theory for ARPES in the case where 
one of the processes can be neglected. We consider the 
case with only double Umklapp scattering in detail and 
present some results for the dimerized case at the end. 
Following we adopt a field theory description of 

the quarter-filled CDW insulator at low energies. The 
small parameter is this approach is the ratio of the gap 
to the bandwidth W. The underlying lattice model may 
be thought of as a quarter-filled Hubbard model with 
additional density-density interactions [^2|Jl^ j9[| . The La- 
grangian density is the sum of two terms describing the 
spin and charge degrees of freedom respectively ("spin- 
charge separation") C = C s + C c , where 

Cs = ^ [V^id^s) 2 +V s (d x $ s ) 2 ] , 

C - = IT-^idr^c) 2 +v c (d x t> s ) 2 ] + Acos(/3$ c ) . (1) 
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Here j3 = \/AK c and the spin sector is gapless |1J] and 
spin rotationally invariant, which fixes the correspond- 
ing Luttinger liquid parameter. The cosine term in the 
charge sector is generated by double Umklapp processes 
and we assume it to be relevant, which implies (3 < 1. We 
will take (3 2 w 0.9 which best fits the estimates given in 
0j for the case of Bechgaard salts. The spin and charge 
velocities v s and v c are parameters of the theory. 

The electronic creation and annihilation operators 
Cj a , Cj jCT (a =t, |) of the underlying lattice model are re- 
lated to the fields appearing in the Lagrangian and their 
corresponding dual fields 



-i f x 
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by 



Cj.o- (x [exTp(ikpx) R a (x) + exp(— ikpx) L a (x)] , (3) 
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where x = jao (do is the lattice spacing) and 



f(|*o+fe c ) ± 5 (<E, s+ e s ) 
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Here ct = ± for up/down spins and r\ a = 77J are Klein 
factors that fulfil {r) ai ri T } = ^a,r- I n order to determine 
the spectral function at zero temperature we evaluate the 
single-electron Green's function 



G { °1(t,x) = (0\R a (T,x)Ri(0)\0). 
We do this by exploiting the factorization ( 

G ( &(T,x) = c (0\Ot(T,x) Oc(0)|0) c 

Ti(* s +e s )(r : x) e ±i(* s +e a )((ii 



(5) 
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where we have introduced the shorthand notation O c for 
the charge part of R a in (||). The subscripts c, s label 
the vacua in the charge and spin sector respectively. The 
spin part of is easily evaluated as the spin sector in 
(0) is a free bosonic theory 



Grr(t,x) = c {0\Ol(r,x) O c (0)|0) c {v. 



(7) 



The large-distance asymptotics of the correlation func- 
tion in the charge sector can be analyzed by going to the 
spectral representation of the corresponding sine-Gordon 
model (SGM). In the repulsive regime (3 > l/v2) the 
only single-particle excitations of the SGM (|l|) are soliton 
and antisoliton. A rough physical picture of what soli- 
tonic excitations look like may be obtained by considering 
a quarter filled extended Hubbard model with sufficiently 
large nearest-neighbour repulsion V |^|| , so that we are 
in the 4fcp CDW insulating phase. One naively obtains 
two ground states as shown in Fig.[l] and the antisoliton 
is then the kink connecting these two ground states. 

• O • O • O • GS1 

O • O • O • O GS2 



• Kink 



FIG. 1. Ground states and antisolitonic kink in a U-V 
extended Hubbard model. 

It is also clear from Fig.[j], the soliton carries an effec- 
tive charge of — e/2. Soliton and antisoliton have massive 
relativistic dispersions, which we as usual parametrize in 
terms of a rapidity variable 

M 

E(9) = M cosh 6 , P(0) = — sinh 9 = m sinh 9 . (8) 

V c 

A basis of scattering states of solitons and antisolitons 
can be constructed by means of the Zamolodchikov- 
Faddeev (ZF) algebra. The ZF algebra is essentially the 



extension of the algebra of creation and annihilation op- 
erators for free fermions or bosons to the case or interact- 
ing particles with factorizable scattering. The ZF algebra 
is based on the knowledge of the exact spectrum and scat- 
tering matrix [ jl5| . For the SGM the ZF operators (and 
their hermitian conjugates) satisfy the following algebra 



Z^(9 1 )Z^(9 2 ) = S^l?(9 1 - 9 2 )Z^(9 2 )Z^(9 1 ) , (9a) 
Zl (9 2 ) = Z\ (9 2 )ZI, {9 1 )Stri (9, - 9 2 ) , (9b) 
Z^(9 1 )ZU9 2 ) = ZU9 2 )S±11(9 2 - 9 l )Z^{9 l ) 



+ (270^(01-02), 



(9c) 



where SJ'J(9) are the known factorizable two-particle 

scattering matrices |1| and ej = — , + for solitons and 
antisolitons respectively. Using the ZF operators a Fock 
space of states can be constructed as follows. The vac- 
uum is defined by 



z £z (9)\n) =0 



(10) 



Multiparticle states are then obtained by acting with 
strings of creation operators Z\ (9) on the vacuum 



|0„ 



l ) en ... £l =Zl(6 n )...Zt i (9 1 )\Q). (11) 



In term of this basis the resolution of the identity is given 
by 



1 = 
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Using this basis of states in the charge sector SGM we 
obtain the following spectral representation 
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The first non- vanishing matrix element (formfactor) in 
( |l2| ) corresponds to the emission of two solitons. Accord- 
ing to Eq.(2.12) of 16 we have 



|<fi|0t(O)|0i, 2 )__| 2 =Z e-^\G{9 1 -9 2 )\ 2 
G{9) = id sinh 0/2 

dt sinh 2 (i[l - z0/tt]) sinh(<[£ - 1]) 
t sinh2t sinh£i cosh 7; 

dt sinh 2 (t/2) sinh(?;[£- 1]) 



x exp 



Ci = exp — 



t sinh 2t sinh £t cosh t 



(13) 



The overall constant Z is calculated in |£6[ and is fixed 
by the following short-distance normalization for r = 

^fv 2 c T 2 + X 2 -> 
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We note that (|lj) implies that Z has dimensions of 

g 2 i 

L 16 2 . In order to restore proper units one needs 
to multiply the correlation function by an appropriate 
power of the lattice spacing. The first subleading contri- 
bution to the correlation function involves three solitons 
and one antisoliton and is expected to be small at low 
energies, so we do not consider it here. The contribu- 
tion of intermediate states involving two solitons to the 
single-electron Green's function is therefore given by 

Z f°° d9- 

T RR\~ 
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Fourier transforming and analytically continuing we find 

dd \G{29)\ 2 



q) = -Z\ 



2vi 



v c q 



where s 2 

v s )/{v c - 



(c(9) + ^c 2 {9) - s 




v 2 q 2 , c(9) = 2Afcosh6>, a = (v c + 



In the limit u. 



v (|T^) simplifies to 



G%hu;,k F + q) 



Z v / v/2n 



\G(26)\ 2 



de- , 

-oo V2Mcosh( 



vq 



1 



2M cosh( 



V(2Mcosh6i) 2 - s 2 
The spectral function for v c — v s is 

Z f b de \G(29)\ 2 y/^ffj 



A RR (uJ,k F + q) = 



vqJo y/s 2 -c 2 (9) 



(17) 



(18) 



where b = arccosh(s/2M) and Z = Z v / 2v/ir 3 . In Fig.| 
we plot Arr(u>, hp + q) as a function of lu for w s = 0.8u c 
and 11 different values of q. As is clear from Fig.|^ the 
spectral function is rather featureless and there are no 
singularities. Furthermore, in contrast to the half-filled 
Mott insulator [ |l9|]T7[ |, there are no dispersing features 
associated separately with with v c and v s . Just above the 
threshold at s 2 — lu 2 — v 2 q 2 = 4M 2 one can approximate 
|G(20)| = |Gi sinh0|, so that for v c = v s 

~~v~ 1 



A RR (u,q) w Z 



(- 



AttM u-vq \ 2M 



(19) 
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FIG. 2. Arr(u, kp+q) as a function of ui/M iovv a — 0.8v c . 
The curves for different q are offset. 

Thus the spectral weight increases linearly with s — 
2M above the threshold. There is a remarkable relation 
between the spectral function and the real part of the 
optical conductivity cr(u)), which for (3 — > 1 and v c = v s 
takes the simple form 



Arr(u, q) 



dx- 



x 5 / 2 a(x) 



const 



— q J 2 m \J lo 2 — v 2 q 2 — x 2 



The tunneling density of states for v c = v s is 



2Z 



L 



arccosh(cL)/2A/ ) 
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2tt v / 2 _ coshl 
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and displays a roughly linear increase after an initial (u> — 
2M) 3 / 2 behaviour just above the threshold at w = 2M. 

Let us now try to relate our results to experiments. 
The ARPES data for PI 23 (Fig. 3 of @) shows a single, 
very broad, dispersing feature that is asymmetric around 
k F and sharpens in the region k > k F . If we interpret 
the underlying increase in intensity in the data as back- 
ground, the signal has a form quite similar to our FigJ^. 
It would be interesting to extract a value Ape for the gap 
from the ARPES data and compare it to gaps seen in op- 
tical measurements A opt and the thermal activation gap 
At extracted e.g. from the dc conductivity. Our theory 
predicts Ape = A opt = 2At- 

Our results may also be of relevance to the Bech- 
gaard salts and in particular to (TMTTF) 2 PF 6 , where 
3D effects appear to be small. However, the ratio of 
the observed gap to the bandwidth is not small, so that 
one should proceed with caution when applying a field- 
theoretical approach. Furthermore, we have neglected 
dimerization, which generates an additional interaction 
// cos in the charge sector of (Q) . The resulting model 
is a double SGM |^0| and not treatable by the method 
used here. If the coupling A of the double Umklapp is 
positive (as is the case for sufficiently strong nearest- 
neighbour repulsion V in the underlying extended Hub- 
bard model) the dimerization does not lead to a qual- 
itative change of the physics described above. If A is 
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negative, the dimerization leads to the confinement of 
the kinks and we expect the physics to be more similar 
to the case without double Umklapp described below. At 
present no definite estimates for U and V for the Bech- 
gaard salts are available. Having said this, let us try 
to compare our results to the experimental findings on 
(TMTTF) 2 PF 6 . The experimental estimates for the gaps 
are A PE = 100 ± 20 meV, A opt w 100 meV and A T « 43 
meV |l8|,fH • This fits well into the small- K c scenario dis- 
cussed above: A opt is twice the single-particle gap be- 
cause the lowest-lying intermediate states that couple to 
the current operator are soliton-antisoliton states plj] . In 
the quarter-filled case, the gap observed in photoemission 
corresponds to the emission of two solitons and therefore 
also occurs at twice the single-particle gap. This is a 
distinct feature of the quarter-filled 4fcp-CDW insulator 
and it looks like there is qualitative agreement with the 
experiments here. 

Let us now turn to the case where we have dimeriza- 
tion, but the interactions are weak so that the double 
Umklapp is irrelevant. The relevant Lagrangian is again 
of the form (Q), but now (3 2 = K c . The main difference to 
the double-Umklapp case is the expression for the elec- 
tron creation operator, which is obtained by replacing 
/3/2 by (3 in (0). The spectral function can then be cal- 
culated by the same methods as before ]l7| , but now 
solitons carry charge e and only a single soliton couples 
to the electron creation operator at low energies. This 
means that now the gap M seen in ARPES is the same 
as the thermal activation gap and half the optical gap. 
In Fig.|| we show the spectral function for v s = 0.Sv c and 
10 different values of q. 
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FIG. 3. Arr(uj, Isf + q) as a function of u>/M with dimer- 
ization but no double Umklapp and v s — 0.8« c . The curves 
for different q are offset and have been smoothed to remove 
the square root threshold singularities. 



For v c q < Q = Mv a / \/v 2 — v 2 The spectral func- 
tion exhibits square root singularities above the thresh- 
old at u> = v/ M 2 + v 2 q 2 . For v c q > Q there are two 
square-root singularities: one above the threshold at 
lo = v s q + Aly/l — (v s /v c ) 2 and a second one when u) 



approaches yj M 2 + v 2 q 2 from below. Their interpreta- 
tion is the same as for the half-filled Mott insulator [jl9| : 
the first singularity correponds to states where the holon 
does not carry momentum, and the second one to states 
where the spinon does not carry momentum. 

This scenario appears to be much less compatible with 
the available ARPES data than the double Umklapp sce- 
nario in the sense that it predicts a rather sharp, dispers- 
ing feature with significant spectral weight around the 
threshold. In the data essentially no dispersing feature 
has been observed, which would appear to be more com- 
patible with the very broad signal shown in Fig. ^| than 
with the one shown in Fig. |^. 
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